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We report on a search for A[J — * p-K~ and A" pK~ (and charge conjugate) decays in pp 
collisions at -^s — 1.96 TeV using 193 pb~^ of data collected by the CDF II experiment at the 
Fermilab Tevatron Collider. Data were collected using a track trigger that has been optimized to 
select tracks belonging to a secondary vertex that is typical of two body charmless decays of b- 
flavored hadrons, including A° baryons. As no Aj signal was observed, we set the upper limits on 
the branching fraction Z3(Ag ph'), where h is if or tt, of 2.3 x 10"^ at 90% C.L. and 2.9 x lO"'^ 
at 95% C.L. 

PACS numbers: 13.30.Eg, 14.20.Mr 



Charmless, hadronic 6-meson decays have been of great 
interest because they provide important information on 
the violation of the combined symmetry operations of 
charge conjugation (C) and parity (P) in the standard 
model of electroweak interactions 0, The first ob- 
servation of charmless hadronic 6-meson decays by the 
CLEO collaboration in 1993 3], and the subsequent real- 
ization that hadronic penguin diagrams dominate some of 
these decays "i^, has since stimulated a substantial body 
of theoretical work In contrast, our present the- 

oretical and experimental knowledge of the correspond- 
ing b baryon decays is rather limited. Measurements of 
branching fractions and CP asymmetries for decays like 



— > pK or pn could provide valuable new insight into 
the hadronic dynamics of 6-hadron decays into charm- 
less final states. In the standard model, the CP-violating 
rate asymmetries in these decays are expected to be large 
compared to the corresponding asymmetries in 6-meson 
decays 

The existence of the A" is well established 0, 0, 
0,0,01, however, no charmless decays have been ob- 
served. We search for A° decaying to p K and p tt. The- 
oretical predictions for their branching ratios lie in the 
range (0.9-1. 2)x 10"*^ forpTr decays and (1.4-1.9) x 10"'^ 
for pK decays jl3| • The current experimental upper limit 
on the branching ratios of these decay modes has been 
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measured by the ALEPH experiment and is 5 x 10~^ at 
90% C.L. [13 ■ The hadronic b trigger of the upgraded 
CoUider Detector at Fermilab (CDF II) selects events 
with track pairs originating from a common displaced 
vertex. A clean signal of charmless hadronic B decays 
has been reconstructed using this trigger . The same 
sample should contain the two-body charmless A^] decays 
in pK and pi:. 

This search uses a 193 ± 12pb~^ data sample 
recorded by the CDF II experiment at the Tevatron pp 
collider with ^/s = 1.96 TeV between February 2002 and 
September 2003. The components of the CDF II detec- 
tor pertinent to this analysis are described briefly below. 
Detailed descriptions can be found elsewhere l21| . Two 
silicon microstrip detectors SVX II [23| and ISL |23| and 
a cylindrical drift chamber COT [23|, immersed in a 1.4 T 
solenoidal magnetic fleld, track charged particles in the 
range \ri\ < 1.0 l^^. The solenoid covers r < 150cm. 
The SVX II provides up to five r — (jj position measure- 
ments, each of ^ 15 fim precision, at radii between 2.5 
and 10.6 cm. The ISL provides one axial and one stereo 
measurement with ~ 20 /im precision, at radii between 
20 and 28 cm, helping to connect the tracks in the COT 
with those in the SVX, and improving the tracking ef- 
ficiency. The COT has 96 measurement layers, between 
40 cm and 137 cm in radius, organized into eight alter- 
nating axial and ±2° stereo superlayers. An additional 
silicon detector, LOO [2^, at radius of 1.3 cm is present 
but is not used in this analysis. 

The events used are selected with a three-level trig- 
ger system. At Level 1, charged tracks in the COT 
transverse plane are reconstructed by a hardware pro- 
cessor (XFT) [231 • The trigger requires two oppositely 
charged tracks with reconstructed transverse momenta 
PT > 2GeV/c and pn + pt2 > 5.5GeV/c. At Level 2, 
the Sihcon Vertex Tracker (SVT) [H associates SVX II 
position measurements with XFT tracks. The impact 
parameter of the track (c?o) with respect to the beam- 
line is measured with 50 /im resolution, which includes 
a ~ 30 /xm contribution from transverse beam size as 
measured in SVT. Requiring two tracks with 100 /xm < 
Mo| < 1.0 mm selects a sample enriched in heavy flavor. 
The two trigger tracks must have an opening angle be- 
tween 20° and 135°. The track pair is also required to 
be consistent with originating from a particle having a 
transverse decay length larger then 200 /im and an im- 
pact parameter less then 140 /im. At Level 3, we fully 
reconstruct the event using the offline software. Candi- 
date trigger tracks are then selected from this improved 
set of tracks by matching them in curvature and (j) to 
tracks reconstructed by the Level 2 trigger. To select 
candidate events, the Level 1 and Level 2 selections are 
then applied to the set of matched tracks, and the invari- 
ant mass, assuming the tracks are pions, is required to 
be within 4 and 7GeV/c'^. 

We normalize the Aj] ph^ branching ratio to 
the branching ratio B {B° ^ Ktt) = (1.85 ±0.11) x 
10~^ I23. The normalization mode has been chosen be- 



cause its decay topology is similar to that of the signal. 
The normalization mode is not well separated from the 
other B h^h'^ decays at CDF, namely B^ tttt, 
S° Kn and S° KK. To obtain the yield of 
Bj^ Ktt we measure the overall B h~^h'~ yield and 
then fit the relative fraction R Ar(S2 -> Ktt)/N{B 
h~^h'~) using an unbinned maximum likelihood fit jl9| . 
The likelihood function has contributions from the signal 
{B h^h'~) and the background. The signal likelihood 
is given by the six distinct decays modes into KK, 
■KTT and K^TT^ . In addition to M^^, the kinematic vari- 
able used is the charged-signed momentum imbalance, 
defined as a = (1 — ^) • qi, where the pi (^2) are the 
modulus of the smaller (larger) momentum of the tracks, 
and qi is the charge sign of the track assigned to pi . 

The relationship between the number of events {N) 
and branching ratios {B) of the signal and normalizing 
mode are given by: 
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where eA (es) is the total efficiency for observing a A° 
(B°) and /a {fd) is the 6-quark hadronization fraction of 
the A° (B°). We use the following values: /a = 0.099 ± 
0.017 and fd = 0.397 ± 0.010 [Ifl. These mean values 
are obtained from measurements at both LEP (see |^ 
[33 ) and CDF [s^, using data samples containing both b 
baryons and mesons and sensitive to pt of the Aj] down 
to 10 GeV/c. The value of the ratio we use is f^lfd = 
0.25 ± 0.04. We estimate the efficiencies using a detailed 
Monte Carlo simulation of the detector and of the trigger 
using GEANT Isl, to generate samples of Ag and B\. 

A blind analysis was performed. The data in the sig- 
nal mass window were hidden and the analysis selections 
optimized without knowledge of their actual impact on 
the result. The background was calculated by fitting the 
invariant mass spectrum and interpolating in the blinded 
signal region. Only after all selection criteria were fixed 
and the systematic uncertainties estimated was the sig- 
nal region unblinded, and the number of events counted 
and compared with the expected background. Potential 
biases in the background estimate, introduced by the cut 
optimization procedure, were avoided by splitting the full 
sample into two statistically independent sub-samples: 
one consisting of even event numbers and the other one 
of odd event numbers. One half of the sample was used 
for the cut optimization described below; the background 
level measured on the other half has been multiplied by 
two to calculate the expected background in the search 
window. 

We select candidate track pairs from the set of offline 
tracks that match trigger tracks based on invariant mass. 
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impact parameter, and transverse decay length of the 
track pair, as well as impact parameter of each track. 
The exact criteria are optimized as discussed below. Fig- 
ure 1 shows the invariant mass distribution after all se- 
lection criteria are applied. The dotted line indicates the 
region that was blinded during the cut optimization. The 
solid line indicates the fit region used to determine the 
expected background level. 




Figure 1: Di pion invariant mass distribution of all the events 
including the search window. The function is the one from 
which we extract the number of B — > h'^h'~ and background 
events. The dashed curve shows the fitted function in the 
part of the mass range that was excluded from the fit. The 
scales of the Monte Carlo distributions of the two signal decay 
modes are arbitrary. The peak in the data is given by the 
B h'^h'~ events. 

We assign the pion mass to all tracks, resulting in slight 
mass shifts between the various 6-hadron decays to two- 
track final states. 

A large Monte Carlo sample including B h^h'^ and 
ph" was used to determine the separation of the 
two mass peaks. For B — > h'^h'" , the mean and rms were 
5258MeV/c2 and 34MeV/c2 . For Ag ph' , the mean 
and rms were 5454MeV/c^ and 60MeV/c^. The separa- 
tion is 196MeV/c^, sufficiently large to make the back- 
ground from B — > h~^h'~ negligible within the search 
window, as can be easily seen from Figure ^ The back- 
ground in the A^ search window is thus predominantly 
combinatoric and can be estimated from the sidebands 
to the right and left of the search window. 

The selection criteria, including the size and position of 
the signal region, were determined from an optimization 
procedure. The sideband regions were defined to include 
those candidates with an invariant mass between 4.800 
and 5.355 GeV/c2 or 5.595 and eCeV/c^. In the opti- 
mization procedure, we take half the candidate events 
in the sideband regions, and maximize a figure of merit 

given by S/ (l.5 + VSK&j 35] where S and BKG rep- 
resent the number of signal and background events, re- 
spectively. The constant in the denominator is chosen 



to favor selections that maximize the sensitivity reach 
at 3(7 significance. This expression reduces to the usual 
SI\/BKG when the background rate is large and to 
5/1.5 when the background is negligible. Hence observ- 
ing that the signal is proportional to the efhciency (ca), 

in the optimization we maximize 



eA/ [X.h + ^JBKG 

where the efficiency has been evaluated using the Monte 
Carlo sample. We simultaneously optimize the cuts on 
the impact parameter of the candidate (^a), its trans- 
verse decay length L^y and the minimum impact param- 
eter of the tracks (min(doi, (io2))- The optimal point 
has been found for |dA| < 50 ^m, Ijxy ^ 400 /im and 
min(dai, do2) > 180 /im. The size and location of the 
mass search window inside the blinded region has been 
optimized according to the same figure of merit and spans 
the mass range between 5.415 GeV/c^ and 5.535 GeV/c^. 

To estimate the expected background, we use the other 
half of the A° sideband sample, and fit it to a sum of a 
Gaussian for the B h^h'~ signal, and various combi- 
nations of exponential and polynomial functions for the 
combinatoric background. The systematic error in the 
yield of i? — > and in the expected combinatoric 

background to the A° — > ph~ signal are estimated from 
the spread of values obtained from different background 
models. Table summarizes these as well as all other 
systematic uncertainties described below. As the central 
value we use the result obtained with the simplest model 
consisting of a Gaussian plus an exponential distribution. 
We arrive at 772 ±31 events for the expected background 
in the A° search window, and 726 ± 82 events for the 
B h^h'^ yield. Uncertainties here include both the 
statistical and systematic errors. 



Affected qty. 


Source 


Syst. Error (%) 


B ^h+h'- yield 


Bkg. shape 


5.7 


Bkg. estimate 


Bkg. shape 


3.3 






3.5 




Window position 


1.2 




Window width 


9 




Lifetime 


3.6 


proton's trigger efficiency 


6 






17 




Overall 


21 


B{Bd ^ Kn) 


5.9 


/A//d 


16 



Table I: List of the relative systematic error contributions to 
the measurement. 

To calculate the 6(Ag -> ph~) we need not only the 
event yields from the data but also the ratio of the effi- 
ciencies es/eA which we evaluate using Monte Carlo sam- 
ples of A" ph~ and B'^ — > Ktt. The efficiency ca was 
evaluated assuming that both A° pn and A° pK 
contribute with the same weight to the signal. We esti- 
mate a systematic uncertainty of 3.5%, allowing for all 
possible values for the ratio of branching fractions. The 
efficiency ratio is also sensitive to the lifetime of the b 



6 



hadron, because the trigger event selection depends on 
the vertex displacement. In the simulation, lifetime val- 
ues from PDG ^ have been used. We varied the lifetime 
values within the experimental uncertainty and observe 
a variation in the efficiency ratio of 3.6%. We quote this 
as a systematic error. 

We assign additional systematic uncertainties due to 
possible discrepancies between Monte Carlo and data 
with regard to invariant mass scale, mass resolution, 
and specific ionization in the COT for different particle 
species. The resulting discrepancies in the mass distri- 
bution are of order a few MeV/c^ and influence signal 
efficiency via the position and width of the search win- 
dow. Varying position and width of the search window to 
reflect the measured differences of data and Monte Carlo 
leads to variations in signal efficiency of 1.2% and 9% 
respectively. 

A third source of systematic error is the variation of 
trigger efficiency with particle species, which arises from 
a different ionization energy loss in the tracking cham- 
ber. We evaluated this effect by adjusting the efficiency 
for pions and kaons using corrections obtained from data. 
As protons and kaons have similar ionization in the mo- 
mentum range of interest, the efficiency for both has been 
corrected in the same way. After the correction, the over- 
all variation of the relative efficiency of 6% was taken as 
the systematic error from this source. 

The main contribution to the systematic error comes 
from the potential difference in pT spectra between b 
mesons and A|]. As the A|j pT spectrum is not well mea- 
sured, we use the h hadron spectrum from and as- 
sume that all hadrons (mesons and baryons) have the 
same spectrum. We compare the efficiency for the inte- 
grated spectrum with the efficiencies for two specific pr 
values. As specific pt values we use the mean of the h 
meson px distribution, and the mean pT of the combi- 
natoric background events below the search window. We 
assign a 17% systematic error based on the spread among 
these three efficiency estimates. 

The value of the efficiency ratio es/eA, corrected for 
the trigger efficiency of different particles, is 1.77 ± 0.37, 
where the error includes both statistical and systematic 
errors. The measured value of the factor A is (3.2±1.0) x 



10^, where statistical and systematic uncertainties are 
included, in addition to uncertainties on B{B^ Ktt) 
and on the production fractions {fd and /A)(Eq. The 
fraction of Ktt {R in eq. ^ is calculated and the 

result is 0.59 ± 0.04. 

The total number of events in the signal region of the 
mass spectrum is 767, consistent within the error with 
the predicted background, 772 ± 31. Because there is no 
excess of signal over the predicted background, we calcu- 
late upper limits on the number of signal events and the 
branching ratio using a Bayesian method with uniform 
prior distribution. This method takes into account the 
effect of statistical and systematic uncertainties. The re- 
sulting upper limits on the number of signal events and 
on B{A^ ph~) are respectively 75 and 2.3 x 10~^ at 
90% C.L and 97 and 2.9 x lO'^ at 95% C.L. This is 
a significant improvement over the previously published 
limit of 5 X 10"^ at 90% C.L for both decay modes fll. 
Substantially more statistics and improved background 
suppression is needed to reach the level of 1 — 2 x 10"^ 
as predicted for the branching fractions in these decays 
in the Standard Model. 
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